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PERIOD THREE IMPLIES CHAOS 


TIEN-YIEN LI AND JAMES A. YORKE 

1. Introduction. The way phenomena or processes evolve or change in time is often described 
by differential equations or difference equations. One of the simplest mathematical situations occurs 
when the phenomenon can be described by a single number as, for example, when the number of 
children susceptible to some disease at the beginning of a school year can be estimated purely as a 
function of the number for the previous year. That is, when the number x n+l at the beginning of the 
n + 1st year (or time period) can be written 

(1.1) x n+l = F(x n ), 

where F maps an interval J into itself. Of course such a model for the year by year progress of the 
disease would be very simplistic and would contain only a shadow of the more complicated 
phenomena. For other phenomena this model might be more accurate. This equation has been used 
successfully to model the distribution of points of impact on a spinning bit for oil well drilling, as 
mentioned in [ 8 , 11 ], knowing this distribution is helpful in predicting uneven wear of the bit. For 
another example, if a population of insects has discrete generations, the size of the n + 1st 
generation will be a function of the nth. A reasonable model would then be a generalized logistic 
equation 

(1.2) x n+i = rx n [l -x nIK]. 

A related model for insect populations was discussed by Utida in [ 10 ]. See also Oster et al [ 14 , 15 ]. 



Fig. 1. For K = 1, r = 3.9, with = .5, the above graph is obtained by iterating Eq. (1.2) 19 times. At right 
the 20 values are repeated in summary. No value occurs twice. While x 2 = .975 and jc, 0 = .973 are close together, 
the behavior is not periodic with period 8 since x w = .222. 


These models are highly simplified, yet even this apparently simple equation (1.2) may have 
surprisingly complicated dynamic behavior. See Figure 1. We approach these equations with the 
viewpoint that irregularities and chaotic oscillations of complicated phenomena may sometimes be 
understood in terms of the simple model, even if that model is not sufficiently sophisticated to allow 
accurate numerical predictions. Lorenz* [ 1 - 4 ] took this point of view in studying turbulent behavior 
in a fascinating series of papers. He showed that a certain complicated fluid flow could be modelled 


985 




986 


TIEN-YIEN LI AND J. A. YORKE 


[December 


by such a sequence x, F(x), F 2 (x), • ■ • , which retained some of the chaotic aspects of the original 
flow. See Figure 2. In this paper we analyze a situation in which the sequence {F n (x)| is non-periodic 
and might be called “chaotic.” Theorem 1 shows that chaotic behavior for (1.1.) will result in any 
situation in which a “population” of size x can grow for two or more successive generations and 
then having reached an unsustainable height, a population bust follows to the level x or below. 



Fig. 2. Lorenz [1] studied the equations for a rotating water-filled vessel which is circularly symmetric about 
its vertical axis. The vessel is heated near the rim and cooled near its center. When the vessel is annular in shape 
and the rotation rate high, waves develop and alter their shape irregularly. From a simplified set of equations 
solved numerically, Lorenz let X„ be in essence the maximum kinetic energy of successive waves. Plotting X n+] 
against X„, and connecting the points, the above graph is obtained. 


In section 3 we give a well-known simple condition which guarantees that a periodic point is 
stable and then in section 4 we quote a result applicable when F is like the one in Figure 2. It implies 
that there is an interval LC / such that for almost every x E /, the set of limit points of the sequence 
(F"(x)} is /oc. 

A number of questions remain unanswered. For example, is the closure of the periodic points an 
interval or at least a finite union of intervals? Other questions are mentioned later. 

Added in proof. May has recently discovered other strong properties of these maps in his 
independent study of how the behavior changes as a parameter is varied [17]. 

2. The main theorem. Let F: /— >/. For x E J, F°(x) denotes x and F n+1 (x) denotes F(F n (x)) 
for n = 0, 1, • • • . We will say p is a periodic point with period n if p E J and p = F n (p) and 
p / F k (p ) for 1 ^ k < n. We say p is periodic or is a periodic point if p is periodic for some n ^ 1 . We 
say q is eventually periodic if for some positive integer m, p = F m (q) is periodic. Since F need not 
be one-to-one, there may be points which are eventually periodic but are not periodic. Our objective 
is to understand the situations in which iterates of a point are very irregular. A special case of our 
main result says that if there is a periodic point with period 3, then for each integer n = 1,2,3, • • • , 
there is a periodic point with period n . Furthermore, there is an uncountable subset of points x in J 
which are not even “asymptotically periodic.” 



1975] 


PERIOD THREE IMPLIES CHAOS 


987 


Theorem 1. Let J be an interval and let F: J J be continuous. Assume there is a point a EJfor 
which the points b = F(a ), c = F 2 (a) and d = F 3 (a), satisfy 

d ^ a < b < c (or d ^ a > b > c). 

Then 

Tl: for every k = 1,2, •• • there is a periodic point in J having period k. 

Furthermore , 

T2: there is an uncountable set S CJ ( containing no periodic points), which satisfies 

the following conditions : 

(A) For evdry p,q E S with p/ q, 

(2.1) limsup \F n (p)-F n (q)\ > 0 

n— »oo 

and 

(2.2) lim inf | F n (p)~ F n (q)j = 0. 

n— *■ °° 

(B) For every p E S and periodic point q EJ, 

limsup \F n (p)~ F n (q)\ > 0. 


Remarks. Notice that if there is a periodic point with period 3, then the hypothesis of the theorem 
will be satisfied. 

An example of a function satisfying the hypotheses of the theorem is F(x) = rx[l - x/K] as in 

(1.2) for r E (3.84, 4] with J = [0, K] and for r > 4, F(x) = max{0, rx[l - x/K]} with J = [0, K]. See 
[2] for a detailed description of iterates of this function for r E [0, 4). The case r = 4 is discussed in 
[6, 7, 12]. 

While the existence of a point of period 3 implies the existence of one of period 5, the converse is 
false. (See Appendix 1). 

We say x EJ is asymptotically periodic if -there is a periodic point p for which 

(2.3) F"(x)-F"(p)^0 as n^ oo. 

It follows from (B) that the set S contains no asymptotically periodic points. We remark that it is 
unknown what the infimum of r is for which the equation (1.2) has points which are not 
asymptotically periodic. 

Proof of Theorem 1. The proof of Tl introduces the main ideas for both Tl and T2. We now give 
the prodf of Tl with necessary lemmas and relegate the tedious proof of T2 to Appendix 2. 

Lemma 0. Let G: I^Rbe continuous, where I is an interval. For any compact interval /, C G(I ) 
there is a compact interval Q Cl such that G(Q) = /,. 

Proof Let /, = [G(p), G(q)], where p,q E I. If p < q, let r be the last point of [p, q] where 
G(r) = G{p) and let s be the l rst point after r where G(s)= G(q). Then G([r, s]) = L. Similar 
reasoning applies when p > q. 

Lemma 1. Let F: J^J be continuous and let {/ n }n=o be a sequence of compact intervals with 
I n C J and 7 n+ ,CF(/ n ) for all n. Then there is a sequence of compact intervals Q n such that 
Q n+ \CQ n Cl 0 and F n (Q n ) = I n for n^ 0. For any x E Q = n Q n we have F"(x) E I n for all n. 
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Proof. Define Q 0 = I 0 . Then F°(Q 0 ) = h- If Qn-i has been defined so that F n_1 (Q n _,) = I n - U then 
I n CF(/ n _,) = F n (Q n -i). By Lemma 0 applied to G = F" on there is a compact interval 
Q n C Q n — i such that F n (Q n )= I n . This completes the induction. 

The technique of studying how certain sequences of sets are mapped into or onto each other is 
often used in studying dynamical systems. For instance, Smale uses this method in his famous 
“horseshoe example” in which he shows how a homeomorphism on the plane can have infinitely 
many periodic points [13]. 

Lemma 2. Let G: J ^ R be continuous. Let I CJ be a compact interval. Assume I C G(I). Then 
there is a point p E I such that G(p) = p. 

Proof. Let I = [(3 0 , j3i]. Choose <*,(/ = 0, 1) in I such that G(a,) = $. It follows a 0 ~ G(a 0 ) = 0 
and a, - G(at) = 0 and so continuity implies G(/3)~ @ must be 0 for some in I. 

Assume d ^ a < b < c as in the theorem. The proof for the case a> b > c is similar and so 
is omitted. Write K = [a, b ] and L = [b,c]. 

Proof of T\: Let k be a positive integer. For k > 1 let {/„} be the sequence of intervals I n = L for 
n = 0,- ",k-2 and I k - X = K, and define I n to be periodic inductively, / n+k = I n torn = 0, 1, 2, • • • . If 
k = 1, let I n = L for all n. 

Let Q n be the sets in the proof of Lemma 1. Then notice that Q k C Q 0 and F k (Q k ) = Qo and so by 
Lemma 2, G = F k has a fixed point p k in Q k . It is clear that p k cannot have period less than k for F; 
otherwise we would need to have F k ~'(p k ) = b , contrary to F k+i (p k ) E L. The point p k is a periodic 
point of period k for F. 

3. Behavior near a periodic point. For some functions F, the asymptotic behavior of iterates of a 
point can be understood simply by studying the periodic points. For 


(3.1) F(x) = ax( 1 - jc) 

a detailed discussion of the points of period 1 and 2 may be found in [1] for a E [0,4] and we now 
summarize some of those results. For a E [0,4], F: [0, 1]— »[0, 1]. 

For a E [0, 1], Jt = 0 is the only point of period 1; in fact, for x E [0, 1], the sequence F n (x )^> 0 as 
n^> oo. 

For a E (1, 3], there are two points of period 1, namely 0 and 1 - a~\ and for x E (0, 1), 
F"(x)^> 1 - a~ x as rc -> <». 

For a > 3 there are also two points of period 2 which we may call p and q and of course F(p) = q 
and F(q) = p. For a E (3, 1 + V6~ 3.449) and x E (0, 1), F 2n (x) converges to either p or q while 
F 2n+1 (x) converges to the other, except for those x for which there is an n for which F"(x) equals 
the point 1 - a~ x of period 1. There are only a countable number of such points so that the behavior 
of {F n (x)| can be understood by studying the periodic points. 

For a > 1 + V6, there are 4 points of period 4 and for a slightly greater than 1 + V6, F 4 "(x) tends 
to one of these 4 unless for some n , F"(x) equals one of the points of period 1 or 2. Therefore we 
may summarize this situation by saying that each point in [0, 1] is asymptotically periodic. 

For those values of a for which each point is asymptotically periodic, it is sufficient to study only 
the periodic points and their “stability properties.” For any function F a point y EJ with period k is 
said to be asymptotically stable if for some interval / = (y - 5, y + 8) we have 

I F k (x)~ y | < |x - y | for all x E I. 

If F is differentiable at the points y, F(y), • • •, F k ~ x (y), there is a simple condition that will guarantee 
this behavior, namely 



< 1 . 
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By the chain rule 

(3.2) F(F-'(y))x-£ F(F‘’(y))x ■ ■ ■ x£ F(y) 

where y n is the nth iterate, F n (y). Therefore y is asymptotically stable if 

flf F(y.) <1, where y, = F‘(y). 

This condition of course guarantees nothing about the limiting behavior of points which do not start 
“near” the periodic point or one of its iterates. The function in Figure 2 which was studied by 
Lorenz has the opposite behavior, namely, where the derivative exists we have 



For such a function every periodic point is “unstable” since for x near a periodic point y of period it, 
the kth iterate F k (x) is further from y than jt is. To see this, approximate F k (x) by 

Fk(y)+ it =y + -£ F ‘<y>u -*]• 

Thus for x near y,\F k (x)-y | is approximately \x - y\\(d/dx)F k (y)\. From (3.2) \(d/dx)F k (y)\ is 
greater than 1. Therefore F n (x) is further from y than x is. 

We do not know when values of a begin to occur for which F in (3.1) has points which are not 
asymptotically periodic. For a = 3.627, F has a periodic point (which is asymptotically stable) of 
period 6 (approx, x = .498). This x is therefore a point of period 3 for F 2 and so Theorem 1 may be 
applied to F 2 . Since F 2 has points which are not asymptotically periodic, the same is true of F. 

In order to contrast the situations in this section with other possible situations discussed in the 
next section, we define the limit set of a point x. The point y is a limit point of a sequence {x„} CJ if 
there is a subsequence x converging to y. The limit set L (x ) is defined to be the set of limit points of 
(F"(x)}. If x is asymptotically periodic, then L(x) is the set {y, F(y), • • •, F^'fy)} for some periodic 
point y of period k. 

4. Statistical properties of (F°(x)}. Theorem 1 establishes the irregularity of the behavior of 
iterates of points. What is also needed is a description of the regular behavior of the sequence 
(F"(x)} when F is piecewise continuously differentiable (as is Lorenz’s function in Figure 2) and 

(4.1) inf > 1 where /, = (x : exists . 

dx { dx 

One approach to describing the asymptotic behavior for such functions is to describe L(x), if 
possible. A second approach, which turns out to be related, is to examine the average behavior of 
{F"(x)}. The fraction of the iterates {x, • • -,F N ~'(x)} of x that are in [a,,a 2 ] will be denoted by 
(t>(x,N,[a l ,a 2 ]). The limiting fraction will be denoted 

<f>(x, [a,, a 2 ]) = lim (t>(x,N,[a„a 2 ]) 

N— >oo 

when the limit exists. The subject of ergodic theory, which studies transformation on general spaces, 
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motivates the following definition. We say g is the density of * (for F) if the limiting fraction satisfies 
(f>(x,[a l ,a 2 ])= \ g(x)dx for all a,,a 2 £/; a,<a 2 . 

J a i 

The techniques for the study of densities use non-elementary techniques of measure theory and 
functional analysis, so that we shall only summarize the results. But their value lies in the fact that 
for certain F almost all jc EJ have the same density. Until recently the existence of such densities 
had not been proved, except for the simplest of functions F. The following result has recently been 
proved: 

Theorem 2. [5]. Let F : /-*/ satisfy the following conditions : 

1) F is continuous. 

2) Except at one point t EJ, F is twice continuously differentiable. 

3) F satisfies (4.1). 

Then there exists a function g : / -> [0, <»), such that for almost all x EJ,g is the density ofx. Also for 
almost all x EJ, L(x) = {y: g(y) > 0} which is an interval. Moreover, the set = {y : g(y) > 0} is an 
interval, and L(x) = /<*, for almost all x. 

The proof makes use of results in [8]. The problem of computationally finding the density is 
solved in [9]. 

A detailed discussion of (3.1) is given in [16], describing how L(x) varies as the parameter a in 

(3.1) varies between 3.0 and 4.0. 

A major question left unsolved is whether (for some nice class of functions F) the existence of a 
stable periodic point implies that almost every point is asymptotically periodic. 

Appendix 1: Period 5 does not imply period 3. In this Appendix we give an example which has a 
fixed point of period 5 but no fixed point of period 3. 

Let F: [1, 5] — > [ 1 , 5], be defined such that F(l) = 3, F( 2) = 5, F(3) = 4, F(4) = 2, F( 5) = 1 and on 
each interval [n, n + 1], 1 ^ n ^4, assume F is linear. Then 

F 3 ([l, 2]) = F 2 ([3, 5]) = F([l, 4]) = [2, 5]. 

Hence, F 3 has no fixed points in [1, 2]. Similarly, F 3 ([2, 3]) = [3, 5] and F 3 ([4, 5]) = [1,4], so neither 
of these intervals contains a fixed point of F 3 . On the other hand, 

F 3 ([3, 4]) = F 2 ([2, 4]) = F([2, 5]) = [1, 5] D [3, 4]. 

Hence, F 3 must have a fixed point in [3, 4]. We shall now demonstrate that the fixed point of F 3 is 
unique and is also a fixed point of F. 

Let p E [3, 4] be a fixed point of F 3 . Then F(p) E [2, 4]. If F(p) E [2, 3], then F 3 (p) would be in 

[1.2] which is impossible since then p could not be a fixed point. Hence F(p)E[3,4] and 
F 2 (p) E [2, 4]. If F 2 (p) E [2, 3] we would have F 3 (p) E [4, 5], an impossibility. Hence p, F{p), F 2 (p) 
are all in [3, 4]. On the interval [3, 4], F is defined linearly and so F(x) = 10 - 2x. It has a fixed point 
10/3 and it is easy to see that F 3 has a unique fixed point, which must be 10/3. Hence there is no point 
of period 3. 

Appendix 2. Proof of T2 of Theorem 1. Let M be the set of sequences M = {M n }n=i of intervals 
with 

(A.l) M n = K or M n CL, and F(M„)DM n+ , 


(A. 2) 


if M n = K then 

n is the square of an integer and M n+U M n +i C L, 
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where K = [a, b] and L = [b, c]. Of course if n is the square of an integer, then n + 1 and n + 2 are 
not, so the last requirement in (A.2) is redundant. For MEM, let P(M, n) denote the number of f s 
in { 1 , • • • , n} for which M, = K. For each r E (3/4, 1) choose M r = {M r „}~=, to be a sequence in M 
such that 

(A.3) iim P(M r , n 2 )ln = r. 

Let M 0 = {M r \ r 6 (3/4, 1)} CM. Then M 0 is uncountable since M r '/M r2 for r, / r 2 . For each 
M r E Mo, by Lemma 1, there exists a point x r with F n (x r ) EM' for all n. Let 5 = {x r : r 6 (3/4, 1)}. 
Then S is also uncountable. For x E S, let P(x, n ) denote the number of /’s in {1, • • •, n} for which 
F‘( x) E K. We can never have F k (x r ) = b, because then would eventually have period 3, contrary 
to (A.2). Consequently P(x r , n ) = P(M r ,n) for all n, and so 

p(x,)= lim P(X n n 2 ) = r 


for all r. We claim that 

(A.4) for p,qES, with p/ q, there exist infinitely many n' s such that F n (p)E K and 
F n {q)EL or vice versa. 

We may assume p(p)>p(q). Then P(p,n)- P(q,n)-^cc, and so there must be infintely many n’s 
such that F n (p)E K and F n (q)EL. 

Since F 2 (b) = d^a and F 2 is continuous, there exists 5 > 0 such that F 2 (x)< ( b + d)l 2 for all 
xE[b-S,b]CK. If pES and F n (p)EK, then (A.2) implies F n+ '(p)EL and F n+2 (p)EL. 
Therefore F"(p)< b - 8. If F n {q)EL, then F n (q)^b so 

I F n (p)~ F n (q)\ > 8. 

By claim (A.4), for any p,q E S, p/ q, it follows 

lim sup | F n (p)~ F n (q) \ g S >0. 

n— *-oc 

Hence (2.1) is proved. This technique may be similarly used to prove (B) is satisfied. 

Proof of 2.2. Since F(b ) = c, F(c) = d ^ a, we may choose intervals [b n , c n ], n = 0, 1,2, • • • , 
such that 

(a) [b, c] = [b°, c°]D [b \ c']D - - - D [b n , c n ]D - - • , 

(b) F(jc) E ( b n ,c n ) for all x E (b n+1 ,c n+l ), 

(c) F{b n+X ) = c n , F(c n+ ') = b n . 

Let A = C\n=o[b n ,c n ], b* = ini A and c* = sup A, then F(b*) = c* and F(c*) = b*, because of 

(c) . 

In order to prove (2.2) we must be more specific in our choice of the sequences M r . In addition to 
our previous requirements on MEM, we will assume that if M k = K for both k = n 2 and ( n + l) 2 
then M k = [b 2n ~ (2i ~ l \b*] for k = n 2 + (2j - 1), M k = [c*,c 2n ~ 2i ] for k = n 2 + 2j where / = 1,- • -,n. 
For the remaining /c’s which are not squares of integers, we assume M k = L. 

It is easy to check that these requirements are consistent with (A.l) and (A.2), and that we can 
still choose M r so as to satisfy (A.3). From the fact that p(x) may be thought of as the limit of the 
fraction of n’s for which F” 2 (jc) E K , it follows that for any r*, r E (3/4, 1) there exist infinitely many 
n such that M k = M r k = K for both k = n 2 and ( n + l) 2 . To show (2.2), let x r E S and x r + E S. Since 
b n b*, c n c* as n for any e > 0 there exists N with | b n - b*\ < e/2, | c n - c*\ < e/2 for 
all n> N. Then, for any n with n > N and M[ = M r k = K for both k = n 2 and ( n + l) 2 , we have 

F n2+l (x r ) E M k = [b 2n ~ x , b*] 
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with k = n 2 + 1 and F n2+l (x r ) and F n2+l (x r *) both belong to [b 2n \b*]. Therefore, | F n2+ \x r )- 
F" 2+1 (jt r *)| < e. Since there are infinitely many n with this property, liminf n ^oo| F n (x r ) - F n (x r *)\ 
= 0 . □ 

Remark. The theorem can be generalized by assuming that F: /-»/? without assuming that 
F(J) CJ and we leave this proof to the reader. Of course F(J) D J would be nonempty since it would 
contain the points a , b , and c, assuming that b , c, and d, are defined. 
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